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[i]  Secondary  circulation,  referring  to  the  motion  relative 
to  a  basic  flow  (geostrophic  and  hydrostatic  balanced), 
occurs  often  in  the  ocean  such  as  deep  convection  and 
circulations  driven  by  fronts  and  eddies.  It  affects  the 
general  circulation  and  the  mass,  heat,  salt,  and  energy 
balance.  The  oceanic  secondary  circulation  is  difficult  to 
measure  directly,  but  is  easy  to  be  identified  by  pseudo- 
vorticity  using  routine  observations.  A  C-vector  method, 
commonly  used  in  atmospheric  mesoscale  moist 
frontogenesis,  is  applied  to  oceanography  for  identifying 
frontal  secondary  circulation  in  Fram  Strait  using 
Conductivity-Temperature-Depth  data  collected  during  a 
large-scale  hydrographic  survey  on  R/V  Valdivia  cruise-54 
of  the  eastern  Greenland  Sea  and  Fram  Strait  from  16  March 
to  5  April  1987.  Possible  application  of  this  method  to 
large-scale  motion  is  also  discussed.  INDEX  TERMS:  4528 
Oceanography:  Physical:  Fronts  and  jets;  4207  Oceanography: 
General:  Arctic  and  Antarctic  oceanography;  4279  Oceanography: 
General:  Upwelling  and  convergences;  4520  Oceanography: 
Physical:  Eddies  and  mesoscale  processes.  Citation:  Chu,  P. 
C.,  C-vector  for  identification  of  oceanic  secondary  circulations 
across  Arctic  Fronts  in  Fram  Strait,  Geophys.  Res.  Lett.,  29(0), 
XXXX,  doi:  10.1 029/2002GL0 15978,  2002. 

1.  Introduction 

[2]  Fram  Strait  (Figure  1)  is  the  only  deep  connection 
between  the  Arctic  Ocean  and  the  rest  of  the  world  ocean 
through  the  Greenland  Sea,  Iceland  Sea,  and  Norwegian  Sea 
(GIN  Sea).  Water  masses  in  Fram  Strait  are  imported  from 
the  neighboring  Atlantic  and  polar  oceans,  and  are  encoun¬ 
tered  in  various  stages  of  modification.  The  North  Atlantic 
Water  (NAW)  is  relatively  warmer  and  saline  (T  >  2°C,  S  > 
34.9  ppt).  The  Polar  Water  (PW)  is  cooler  and  fresher  (T< 
0°C,  S  <  34.7  ppt)  [van  Aken  et  al.,  1991].  Analyzing  1984 
marginal  ice  zone  experiment  observations  in  Fram  Strait, 
Gas  card  et  al.  [1988]  found  that  the  West  Spitsbergen 
Current  and  the  East  Greenland  Current  are  two  main 
generators  for  eddies  in  Fram  Strait. 

[3]  In  the  transition,  different  water  masses  interface  and 
form  frontal  zones  that  not  only  separate  water  bodies  with 
different  hydrographic  characteristics  but  also  the  regional 
biological  systems.  The  different  water  masses  encountered 
in  the  GIN  Sea  and  Fram  Strait  often  form  fronts.  The  Arctic 
front  or  frontal  zone  is  oriented  more  or  less  meridionally. 
It  separates  the  wann  and  salty  NAW  in  the  Norwegian  and 
West  Spitsbergen  currents  from  the  cooler  and  fresher  Arctic 
Water  (AW)  [ Dietrich ,  1969].  After  analyzing  hydrographic 


data  along  74°45'N  in  February  1989  during  cruise  VA78  of 
R/V  Valdivia,  van  Aken  et  al.  [1991]  identified  four  Arctic 
fronts  south  of  Fram  Strait.The  physical-biological  effect 
arises  from  the  significant  vertical  component  of  the  3D 
ageostrophic  flow  (or  called  the  secondary  circulation) 
associated  with  the  fronts,  where  vertical  upward  motion 
may  act  as  a  fertilizer  of  the  upper  water  column.  A  better 
knowledge  of  the  secondary  circulation  is  therefore  crucial. 
However,  direct  measurement  of  the  secondary  circulation 
in  the  ocean  is  difficult. 

[4]  Use  of  vorticity  for  depicting  3D  secondary  circula¬ 
tion  is  popular  in  meteorology.  For  example,  vertical  vor¬ 
ticity  component  can  be  used  to  describe  gyres  and  eddies 
and  horizontal  vorticity  component  can  be  used  to  depict  the 
vertical  circulation.  The  Q-vector  concept  proposed  by 
Hoskins  et  al.  [1978]  and  related  analysis  methods  on  the 
base  of  quasi-geostrophic  theory  have  been  very  useful  in 
understanding  and  diagnosing  the  frontal  secondary  circu¬ 
lation.  However,  the  barotropic  part  of  the  rotational  ageo¬ 
strophic  flow  is  excluded  in  the  Q-vector  equation.  To 
overcome  this  weakness,  Xu  [1992]  proposed  a  C-vector 
to  explain  atmospheric  mesoscale  moist  frontogenesis. 
However,  the  C-vector  analysis  hasn’t  been  applied  to 
physical  oceanography  yet.  This  study  shows  the  usefulness 
of  the  C-vector  in  diagnosing  three-dimensional  oceanic 
secondary  circulations  from  observational  data. 


2.  C-Vector 

[5]  Let  (x,  y)  be  the  horizontal  coordinates  and  z  the 
vertical  coordinate.  With  the  Boussinesq  approximation,  the 
geostrophic  velocity  (wg,  vg)  is  related  to  the  density  through 
the  thermal  wind  relation, 
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where  p  =  p  -  p0,  and  (p,  p0)  are  the  in-situ  density  and 
characteristic  density.  The  variable  b  is  usually  called  the 
buoyancy.  The  total  flow,  V  =  (u,  v,  w),  is  decomposed  into 
geostrophic  and  ageostrophic  parts:  V  =  Vg  +  Vag.  If  the 
advection  of  momentum  and  buoyancy  is  dominated  by  the 
geostrophic  advection  (i.e.,  quasi-geostrophic  system), 


V.VVRiVg.  VVg,  V  •  Vb  «  Vg  •  Vb,  (2) 


the  ageostrophic  velocity  is  determined  by 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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Longitude 


(c)  buoyancy  flux  (B).  In  the  upper  ocean,  the  last  two 
factors  are  mainly  caused  by  the  surface  wind  stress  and 
buoyancy  flux.  As  pointed  by  Xu  [1992],  the  C-vector  is  the 
ageostrophic  vortex  line  (Figure  2)  whose  horizontal 
components  ( Cx ,  Cy)  represent  the  secondary  circulation  in 
vertical  cross-section. 

[7]  Observations  show  generally  well-mixed  upper  ocean 
by  turbulent  motion.  The  transition  between  the  turbulent 
mixed  layer  and  the  stratified  water  below  is  a  thin  entrain¬ 
ment  zone  with  large  gradients  of  density  and  velocity.  All 
turbulent  fluxes  are  usually  assumed  to  vanish  below  the 
ocean  mixed  layer.  The  mixed  layer  depth  h  is  defined  as  the 
depth  above  which  temperature,  salinity,  or  velocity  (geo- 
strophic  plus  ageostrophic)  is  vertically  uniform  to  certain 
critical  value.  The  vertical  uniformity  leads  to  the  bulk 
model  parameterization  [e.g.,  Price  et  al.,  1986;  Chu  et 
al.,  1990;  Chu  and  Garwood,  1991]  of  turbulent  fluxes  ( X , 
Y,B ) 


Figure  1.  Geography  of  East  Fram  Strait,  coordinate 

system,  and  CTD  stations  of  R/V  Valdivia  cruise  54  in  Fram  (V,  Y)  =  ( tx ,  ry)  +  |  (rx.  ry)  -  (X,  Y)  h  "1 , 

Strait  from  March  16  to  April  5,  1987. 

B  =  B0  +  (B0  -  5_*)  j,  for  z  >  -h,  (8a) 
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where  ( X ,  Y)  and  B  are  the  vertical  turbulent  momentum  and 
buoyancy  fluxes  (downward  positive),  and  N  is  the  Brunt- 
Vaisala  frequency. 


[6]  Cross  derivatives  among  (3)-(5)  lead  to  the  definition 
of  pseudovorticity  [Xu,  1992], 

|(^) 

(6a) 

-lx(N^a)=2Cy, 

(6b) 

-^-y{f2ua)=2C„ 

(6c) 

(X,  7,5)  ~  0  for  z  <  —h,  (8b) 

where  (tx,  Ty)  is  the  surface  wind  stress,  and  B0  is  the 
surface  buoyancy  flux.  ( X ,  Y)_h  and  B_h  are  turbulent  fluxes 
at  the  base  of  mixed  layer,  which  are  computed  from  the 
surface  fluxes.  From  (8a)  and  (8b)  we  obtain 


d2X 

dz2 


=  0, 


(9) 


for  the  whole  water  column.  In  computing  the  horizontal 
components  of  the  C-vector  ( Cx ,  Cv),  the  turbulent 
momentum  flux  can  be  neglected  for  the  whole  water 
column,  and  the  turbulent  buoyancy  flux  can  be  neglected 
for  the  water  column  below  the  surface  mixed  layer.  Note 
that  non-unifonn  currents  may  exist  in  the  upper  ocean  such 
as  the  Ekman  layer,  where  (9)  does  not  hold.  Flowever, 
equation  (9)  may  be  used  below  the  Ekman  layer. 

[8]  Since  B  is  a  linear  function  of  z  in  the  surface  mixed 
layer  [see  (8a)],  the  contribution  of  B  to  (Cx,  Cy )  in  the 
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are  components  of  C-vector.  Thus,  the  pseudovorticity  of 
the  secondary  circulation  is  determined  by  three  forcing 
factors:  (a)  geostrophic  forcing  (i.e.,  distinct  water  masses 
acrossthefront),— fd(ug,v^ld{y,z),— fd(ug,vg)/d(z,x),— fd(ug, 
—fd(ug,  vg)/<9(x,  y),  (b)  turbulent  momentum  flux  ( X ,  Y),  and 


CHU:  C-VECTOR  FOR  IDENTIFICATION  OF  OCEANIC  SECONDARY  CIRCULATIONS 


X  -  3 


STATION 


Figure  3.  Potential  density  excess  referred  to  500  db  along  the  north  cross-section  in  Fram  Strait  (from  Stations  44-51). 


surface  mixed  layer  is  depth-independent  [see  (8a)].  Since 
B_h  is  computed  from  the  surface  turbulent  fluxes  (tx,  tv) 
and  Bq  (bulk  models),  the  atmosphere  may  control  the  scale 
of  the  horizontal  variability  of  B_h  and  B0.  The  atmosphere 
has  most  of  its  energy  in  scales  of  several  hundreds  of 
kilometers  and  longer,  it  is  reasonable  to  assume  that 


dz  \i9y ) 


=  0, 


d_  fdB\ 
Fz\dx) 


=  o, 


(10) 


across  the  Arctic  front,  and  therefore  the  horizontal  C-vector 
components  are  computed  from  the  geostrophic  forcing 
only.  The  geostrophic  forcing  has  been  recognized  as  the 
major  factor  to  cause  the  atmospheric  frontal  secondary 
circulation  [. Hoskins  et  al.,  1978;  Xu,  1992],  It  may  also  be 
important  for  oceanic  frontal  secondary  circulation. 


3.  Data 

[9]  Conductivity-Temperature-Depth  (CTD)  data,  col¬ 
lected  during  a  large-scale  hydrographic  survey  on  RV/ 
VALDIVIA  cruise-54  of  the  eastern  Greenland  Sea  and 
Fram  Strait  from  16  March  to  5  April  1987  [Quadfasel  and 
Ungewiss,  1988],  are  used  to  illustrate  the  advantage  of 
using  the  C-vector  method  in  analyzing  oceanic  secondary 
circulation.  The  major  task  of  this  cruise  was  to  map  the 
vertical  distribution  of  temperature,  salinity,  and  dissolved 
oxygen  in  the  Greenland  Sea  as  a  measure  of  the  large-scale 
circulation  and  transport.  A  secondary  objective  was  to 
search  active  convection  events.  Along  seven  sections  a 
total  of  73  CTD  profiles  were  taken  (Figure  1).  Four  of 
these  sections  crossed  the  Arctic  Front  that  separates  the 
Greenland  Sea  gyres  from  the  warm  and  saline  northward 
flowing  Norwegian  Atlantic  Current  and  the  West  Spitsber¬ 
gen  Current.  The  sections  were  designed  to  form  three 
closed  boxes  to  allow  calculation  of  transport  budgets. 
Usual  station  spacing  was  56  km  except  along  Fram  Strait 
section  at  78°  25'N  and  across  the  Hovgaard  Fracture  Zone, 


where  the  spacing  was  decreased  to  less  than  28  km.  All 
CTD  profiles  were  run  to  within  5  m  of  the  bottom. 

4.  Potential  Density 

[10]  The  potential  density  excess  referred  to  500  db 
computed  from  the  CTD  data  shows  the  existence  of 
multi-frontal  zones  in  Fram  Strait.  For  example,  three  Arctic 
fronts  are  identified  from  the  north  cross  section  (Stations 
44  to  51)  (Figure  3):  (a)  eastern  front  occurring  near  the 
west  coast  of  Spitsbergen  (Stations  49-51),  (b)  shallow 
western  front  (above  160  m  depth)  occurring  near  0° 
longitude  (Stations  44-46)  in  the  Hovgaard  Fracture  Zone, 
and  (c)  a  deep  mid-front  (Stations  47-48)  in  the  north 
Knipovich  Ridge.  In  the  south  cross-section  in  Fram  Strait 
(Stations  26-41,  Figure  4),  only  two  fronts  can  be  identified 
with  a  strong  and  shallow  front  in  the  Boreas  Basin 
(Stations  34-36)  and  a  weak  deep  front  near  the  west  coast 
of  Spitsbergen  (Stations  26-30). 

[11]  Usually,  upward  (downward)  bending  of  isopycnals 
is  used  to  identify  upward  (downward)  motion.  For  exam¬ 
ple,  in  Figure  3  the  isopycnals  bend  upward  near  Stations 
46-47,  49,  and  downward  near  Stations  45  (to  250  m 
depth),  48  (surface  to  1000  m  depth),  and  51.  Since  such 
identification  is  very  qualitative,  it  is  hard  to  sketch  the 
secondary  circulation  pattern  and  strength.  Is  it  possible  to 
get  quantitative  information  using  the  same  density  data? 
The  answer  is  ’yes’  because  the  secondary  circulation  can 
be  represented  by  the  horizontal  component  of  the  C-vector 
(i.e.,  horizontal  pseudovorticity). 

5.  Horizontal  Pseudovorticity 

[12]  Taking  2500  db  as  the  reference  level,  the  geo¬ 
strophic  current  (ug,  vg)  is  computed  from  density,  and  then 
the  horizontal  components  of  the  C-vector  are  computed 
from  (Mg,  vg).  Figure  5  shows  the  x-component  of  the  non- 
dimensional  pseudovorticity  ( CJf 2)  along  the  north  cross 
section  (Stations  44  to  51).  Looking  toward  north  the 
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Figure  4.  Potential  density  excess  referred  to  500  db  along  the  south  cross-section  in  Fram  Strait  (from  Stations  26-41). 


positive  (negative)  values  of  CJf2  imply  a  clockwise  (anti¬ 
clockwise)  circulation.  Two  clockwise  and  two  anticlock¬ 
wise  secondary  circulations  are  identified.  Among  them,  the 
clockwise  secondary  circulations  are  deeper  than  the  anti¬ 
clockwise  secondary  circulations  (depth  less  than  60  m). 
Note  that  in  this  study,  the  geostrophic  current  is  computed 
as  2500  db  to  be  assumed  as  the  level  of  no  motion.  This 
may  affect  the  C-vector  computation  especially  in  large 
horizontally  sheared  barotropic  current  such  as  in  Fram 
Strait  [Fahrbach  et  al.,  2001],  However,  Figure  5  keeps 
almost  the  same  as  2000  db  is  chosen  as  the  reference  level. 

[13]  The  most  striking  feature  is  the  existence  of  a 
shallow  (surface  to  60  m),  strong  anticlockwise  secondary 
circulation  with  a  minimum  value  of  CJf 2  around  —9.01, 


located  near  the  Knipovich  Ridge  (Stations  46-48).  We 
may  call  it  the  Knipovich  cell.  Its  upward  and  downward 
branches  are  connected  to  two  clockwise  secondary  cells 
from  the  east  (upward  branch)  and  west  (downward  branch). 
Below  the  Knipovich  cell  between  340  and  1000  m  depth, 
there  is  a  very  weak  clockwise  vorticity  (a  maximum  value 
of  CJf 2  around  0.01). 

[14]  East  of  the  Knipovich  cell,  a  deep  (surface  to  1000 
m),  clockwise  secondary  circulation  with  a  maximum 
pseudovorticity  of  1 .95  is  identified  near  the  West  Spitsber¬ 
gen  slope  (Stations  48-50).  We  may  call  it  the  West 
Spitsbergen  cell.  Its  downward  branch  follows  the  slope. 
Its  upward  branch  connects  to  the  Knipovich  cell.  West  of 
the  Knipovich  cell,  a  relatively  shallow  (to  250  m  depth), 


STATION 

Figure  5.  Nondimensional  horizontal  pseudovorticity,  CJf2,  along  the  north  cross-section  in  Fram  Strait  (from  Stations 
44-51). 
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clockwise  secondary  circulation  with  a  maximum  pseudo- 
vorticity  of  3.16  is  identified  in  the  Hovgaard  Fracture  Zone 
(Stations  44-46).  Downward  branch  connects  to  the  Kni- 
povich  cell  (Stations  46-48).  The  upward  branch  is  located 
at  Stations  44-45.  We  may  call  it  the  Hovgaard  cell.  Since 
the  horizontal  pseudovorticity  ( CJf 2)  represents  the  secon¬ 
dary  circulation  around  the  x-axis.  The  larger  the  CJf2,  the 
stronger  the  secondary  circulation  is.  Strong  upwelling  is 
identified  between  the  Knipovich  Ridge  and  the  West 
Spitsbergen  slope. 


6.  Possible  Application  to  Large-Scale  Motion 

[15]  For  large-scale  motion,  the  horizontal  variability  of 
B_h  and  B0  cannot  be  assumed  as  zero,  and  (10)  is  not  valid. 
However,  Equation  (9)  is  still  valid  since  it  does  not  depend 
on  the  horizontal  scale.  The  C-vector  for  the  large-scale 
motion  can  be  written  by 


r  f9{«g,vg)  1  d  fdB\ 

1  d(y,  z)  2dz\dy)' 

r  .d(»g,vg)  1  d  fdB\ 

y  1  d{z,x)  2  dz\dx)' 

_  d(ug,vg)  f  a  fax  ar\  por 
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(11a) 

(lib) 

(11c) 


where  the  thermohaline  circulation  can  be  represented  by 
the  horizontal  components  (Cx,  Cy)  and  the  wind-driven 
circulation  can  be  represented  by  the  vertical  component  Cz. 
The  overturning  streamfunction  can  be  obtained  from  (Cx, 
Cy)  by  solving  the  Poisson  equation  with  the  C-vector 
components  as  the  forcing  terms. 


7.  Conclusions 

[16]  This  paper  illustrates  the  usefulness  of  the  C-vector 
method  in  identifying  the  secondary  circulations  in  Fram 
Strait  frontal  zone.  However,  this  method  might  be  used  to 
diagnose  the  large-scale  motion  such  as  the  thermohaline 


circulation  represented  by  the  horizontal  components  (Cx, 
Cy)  and  the  wind-driven  circulation  represented  by  the 
vertical  component  Cz.  The  overturning  streamfunction 
can  be  obtained  from  (Cx,  Cy)  by  solving  the  Poisson 
equation  with  the  C-vector  components  as  the  forcing  terms. 
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